Genomic insights into Spiroplasma endosymbionts that induce male-killing and protective phenotypes in the pea aphid.

The endosymbiotic bacteria Spiroplasma (Mollicutes) infect diverse plants and arthropods, and some of which induce male killing, where male hosts are killed during development. Male-killing Spiroplasma strains belong to either the phylogenetically distant Citri-Poulsonii or Ixodetis groups. In Drosophila flies, Spiroplasma poulsonii induces male killing via the Spaid toxin. While Spiroplasma ixodetis infects a wide range of insects and arachnids, little is known about the genetic basis of S. ixodetis-induced male killing. Here, we analyzed the genome of S. ixodetis strains in the pea aphid Acyrthosiphon pisum (Aphididae, Hemiptera). Genome sequencing constructed a complete genome of a male-killing strain, sAp269, consisting of a 1.5 Mb circular chromosome and an 80 Kb plasmid. sAp269 encoded putative virulence factors containing either ankyrin repeat, ovarian tumor-like deubiquitinase, or ribosome inactivating protein domains, but lacked the Spaid toxin. Further comparative genomics of Spiroplasma strains in A. pisum biotypes adapted to different host plants revealed their phylogenetic associations and the diversity of putative virulence factors. Although the mechanisms of S. ixodetis-induced male killing in pea aphids remain elusive, this study underlines the dynamic genome evolution of S. ixodetis and proposes independent acquisition events of male-killing mechanisms in insects.

Identification of a viral gene essential for the genome replication of a domesticated endogenous virus in ichneumonid parasitoid wasps.

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposotor didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis.

Chromosome-Level Assembly and Annotation of the Pearly Heath Coenonympha arcania Butterfly Genome.

We present the first chromosome-level genome assembly and annotation of the pearly heath Coenonympha arcania, generated with a PacBio HiFi sequencing approach and complemented with Hi-C data. We additionally compare synteny, gene, and repeat content between C. arcania and other Lepidopteran genomes. This reference genome will enable future population genomics studies with Coenonympha butterflies, a species-rich genus that encompasses some of the most highly endangered butterfly taxa in Europe.

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid.

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

MTG-Link: leveraging barcode information from linked-reads to assemble specific loci.

BACKGROUND: Local assembly with short and long reads has proven to be very useful in many applications: reconstruction of the sequence of a locus of interest, gap-filling in draft assemblies, as well as alternative allele reconstruction of large Structural Variants. Whereas linked-read technologies have a great potential to assemble specific loci as they provide long-range information while maintaining the power and accuracy of short-read sequencing, there is a lack of local assembly tools for linked-read data. RESULTS: We present MTG-Link, a novel local assembly tool dedicated to linked-reads. The originality of the method lies in its read subsampling step which takes advantage of the barcode information contained in linked-reads mapped in flanking regions. We validated our approach on several datasets from different linked-read technologies. We show that MTG-Link is able to assemble successfully large sequences, up to dozens of Kb. We also demonstrate that the read subsampling step of MTG-Link considerably improves the local assembly of specific loci compared to other existing short-read local assembly tools. Furthermore, MTG-Link was able to fully characterize large insertion variants and deletion breakpoints in a human genome and to reconstruct dark regions in clinically-relevant human genes. It also improved the contiguity of a 1.3 Mb locus of biological interest in several individual genomes of the mimetic butterfly Heliconius numata. CONCLUSIONS: MTG-Link is an efficient local assembly tool designed for different linked-read sequencing technologies. MTG-Link source code is available at https://github.com/anne-gcd/MTG-Link and as a Bioconda package.

Genome-wide identification of lncRNAs associated with viral infection in Spodoptera frugiperda.

The role of lncRNAs in immune defence has been demonstrated in many multicellular and unicellular organisms. However, investigation of the identification and characterization of long non-coding RNAs (lncRNAs) involved in the insect immune response is still limited. In this study, we used RNA sequencing (RNA-seq) to investigate the expression profiles of lncRNAs and mRNAs in the fall armyworm Spodoptera frugiperda in response to virus infection. To assess the tissue- and virus-specificity of lncRNAs, we analysed and compared their expression profiles in haemocytes and fat body of larvae infected with two entomopathogenic viruses with different lifestyles, i.e. the polydnavirus HdIV (Hyposoter didymator IchnoVirus) and the densovirus JcDV (Junonia coenia densovirus). We identified 1883 candidate lncRNAs, of which 529 showed differential expression following viral infection. Expression profiles differed considerably between samples, indicating that many differentially expressed (DE) lncRNAs showed virus- and tissue-specific expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and target prediction analyses indicated that DE-LncRNAs were mainly enriched in metabolic process, DNA replication and repair, immune response, metabolism of insect hormone and cell adhesion. In addition, we identified three DE-lncRNAs potentially acting as microRNA host genes, suggesting that they participate in gene regulation by producing miRNAs in response to virus infection. This study provides a catalogue of lncRNAs expressed in two important immune tissues and potential insight into their roles in the antiviral defence in S. frugiperda. The results may help future in-depth functional studies to better understand the biological function of lncRNAs in interaction between viruses and the fall armyworm.

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies.

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

The evolutionary process of invasion in the fall armyworm (Spodoptera frugiperda).

The fall armyworm (FAW; Spodoptera frugiperda) is one of the major agricultural pest insects. FAW is native to the Americas, and its invasion was first reported in West Africa in 2016. Then it quickly spread through Africa, Asia, and Oceania, becoming one of the main threats to corn production. We analyzed whole genome sequences of 177 FAW individuals from 12 locations on four continents to infer evolutionary processes of invasion. Principal component analysis from the TPI gene and whole genome sequences shows that invasive FAW populations originated from the corn strain. Ancestry coefficient and phylogenetic analyses from the nuclear genome indicate that invasive populations are derived from a single ancestry, distinct from native populations, while the mitochondrial phylogenetic tree supports the hypothesis of multiple introductions. Adaptive evolution specific to invasive populations was observed in detoxification, chemosensory, and digestion genes. We concluded that extant invasive FAW populations originated from the corn strain with potential contributions of adaptive evolution.

Host-plant adaptation as a driver of incipient speciation in the fall armyworm (Spodoptera frugiperda).

BACKGROUND: Divergent selection on host-plants is one of the main evolutionary forces driving ecological speciation in phytophagous insects. The ecological speciation might be challenging in the presence of gene flow and assortative mating because the direction of divergence is not necessarily the same between ecological selection (through host-plant adaptation) and assortative mating. The fall armyworm (FAW), a major lepidopteran pest species, is composed of two sympatric strains, corn and rice strains, named after two of their preferred host-plants. These two strains have been hypothesized to undergo incipient speciation, based on (i) several lines of evidence encompassing both pre- and post-zygotic reproductive isolation, and (ii) the presence of a substantial level of genetic differentiation. Even though the status of these two strains has been established a long time ago, it is still yet to be found whether these two strains indeed exhibit a marked level of genetic differentiation from a large number of genomic loci. Here, we analyzed whole genome sequences from 56 FAW individuals either collected from pasture grasses (a part of the favored host range of the rice strain) or corn to assess the role of host-plant adaptation in incipient speciation. RESULTS: Principal component analysis of whole genome data shows that the pattern of divergence in the fall armyworm is predominantly explained by the genetic differentiation associated with host-plants. The level of genetic differentiation between corn and rice strains is particularly marked in the Z chromosome. We identified one autosomal locus and two Z chromosome loci targeted by selective sweeps specific to rice strain and corn strain, respectively. The autosomal locus has both increased DXY and FST while the Z chromosome loci had decreased DXY and increased FST. CONCLUSION: These results show that the FAW population structure is dominated by the genetic differentiation between corn and rice strains. This differentiation involves divergent selection targeting at least three loci, which include a locus potentially causing reproductive isolation. Taken together, these results suggest the evolutionary scenario that host-plant speciation is a driver of incipient speciation in the fall armyworm.

Spodoptera littoralis genome mining brings insights on the dynamic of expansion of gustatory receptors in polyphagous noctuidae.

The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species.

Multi-Omic Investigation of Low-Nitrogen Conditional Resistance to Clubroot Reveals Brassica napus Genes Involved in Nitrate Assimilation.

Nitrogen fertilization has been reported to influence the development of clubroot, a root disease of Brassicaceae species, caused by the obligate protist Plasmodiophora brassicae. Our previous works highlighted that low-nitrogen fertilization induced a strong reduction of clubroot symptoms in some oilseed rape genotypes. To further understand the underlying mechanisms, the response to P. brassicae infection was investigated in two genotypes "Yudal" and HD018 harboring sharply contrasted nitrogen-driven modulation of resistance toward P. brassicae. Targeted hormone and metabolic profiling, as well as RNA-seq analysis, were performed in inoculated and non-inoculated roots at 14 and 27 days post-inoculation, under high and low-nitrogen conditions. Clubroot infection triggered a large increase of SA concentration and an induction of the SA gene markers expression whatever the genotype and nitrogen conditions. Overall, metabolic profiles suggested that N-driven induction of resistance was independent of SA signaling, soluble carbohydrate and amino acid concentrations. Low-nitrogen-driven resistance in "Yudal" was associated with the transcriptional regulation of a small set of genes, among which the induction of NRT2- and NR-encoding genes. Altogether, our results indicate a possible role of nitrate transporters and auxin signaling in the crosstalk between plant nutrition and partial resistance to pathogens.

RNA sequencing data for responses to drought stress and/or clubroot infection in developing seeds of Brassica napus.

Oilseed rape (Brassica napus L.) is the third largest oil crop worldwide. Like other crops, oilseed rape faces unfavorable environmental conditions resulting from multiple and combined actions of abiotic and biotic constraints that occur throughout the growing season. In particular drought severely reduces seed yield but also impacts seed quality in oilseed rape. In addition, clubroot disease, caused by the pathogen Plasmodiophora brassicae, limits the yield of the oilseed rape crops grown in infected areas. Clubroot induces swellings or galls on the roots that decrease the flow of water and nutrients within the plant. Furthermore, combinations of different stresses lead to complex plant responses that can not be predicted by the simple addition of individual stress responses. Indeed, an abiotic constraint can either reduce or stimulate the plant response to a pathogen or pest. Transcriptome datasets from different conditions are key resources to improve our knowledge of environmental stress-resistance mechanisms in plant organs. Here, we describe a RNA-seq dataset consisting of 72 samples of immature B. napus seeds from plants grown either under drought, infected with P. brassicae, or a combination of both stresses. A total of 67.6 Gb of transcriptome paired-end reads were filtered, mapped onto the B. napus reference genome Darmor-bzh and used for identification of differentially expressed genes and gene ontology enrichment. The raw reads are available under accession PRJNA738318 at NCBI Sequence Read Archive (SRA) repository. The dataset is a resource for the scientific community exploring seed plasticity.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus.

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

Proteo-Trancriptomic Analyses Reveal a Large Expansion of Metalloprotease-Like Proteins in Atypical Venom Vesicles of the Wasp Meteorus pulchricornis (Braconidae).

Meteorus pulchricornis (Ichneumonoidea, Braconidae) is an endoparasitoid wasp of lepidopteran caterpillars. Its parasitic success relies on vesicles (named M. pulchricornis Virus-Like Particles or MpVLPs) that are synthesized in the venom gland and injected into the parasitoid host along with the venom during oviposition. In order to define the content and understand the biogenesis of these atypical vesicles, we performed a transcriptome analysis of the venom gland and a proteomic analysis of the venom and purified MpVLPs. About half of the MpVLPs and soluble venom proteins identified were unknown and no similarity with any known viral sequence was found. However, MpVLPs contained a large number of proteins labelled as metalloproteinases while the most abundant protein family in the soluble venom was that of proteins containing the Domain of Unknown Function DUF-4803. The high number of these proteins identified suggests that a large expansion of these two protein families occurred in M. pulchricornis. Therefore, although the exact mechanism of MpVLPs formation remains to be elucidated, these vesicles appear to be "metalloproteinase bombs" that may have several physiological roles in the host including modifying the functions of its immune cells. The role of DUF4803 proteins, also present in the venom of other braconids, remains to be clarified.

Plant-phenotypic changes induced by parasitoid ichnoviruses enhance the performance of both unparasitized and parasitized caterpillars.

There is increasing awareness that interactions between plants and insects can be mediated by microbial symbionts. Nonetheless, evidence showing that symbionts associated with organisms beyond the second trophic level affect plant-insect interactions are restricted to a few cases belonging to parasitoid-associated bracoviruses. Insect parasitoids harbour a wide array of symbionts which, like bracoviruses, can be injected into their herbivorous hosts to manipulate their physiology and behaviour. Yet, the function of these symbionts in plant-based trophic webs remains largely overlooked. Here, we provide the first evidence of a parasitoid-associated symbiont belonging to the group of ichnoviruses which affects the strength of plant-insect interactions. A comparative proteomic analysis shows that, upon parasitoid injection of calyx fluid containing ichnovirus particles, the composition of salivary glands of caterpillars changes both qualitatively (presence of two viral-encoded proteins) and quantitatively (abundance of several caterpillar-resident enzymes, including elicitors such as glucose oxidase). In turn, plant phenotypic changes triggered by the altered composition of caterpillar oral secretions affect the performance of herbivores. Ichnovirus manipulation of plant responses to herbivory leads to benefits for their parasitoid partners in terms of reduced developmental time within the parasitized caterpillar. Interestingly, plant-mediated ichnovirus-induced effects also enhance the performances of unparasitized herbivores which in natural conditions may feed alongside parasitized ones. We discuss these findings in the context of ecological costs imposed to the plant by the viral symbiont of the parasitoid. Our results provide intriguing novel findings about the role played by carnivore-associated symbionts on plant-insect-parasitoid systems and underline the importance of placing mutualistic associations in an ecological perspective.

Global patterns in genomic diversity underpinning the evolution of insecticide resistance in the aphid crop pest Myzus persicae.

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

Transcriptome Profiling of Starvation in the Peripheral Chemosensory Organs of the Crop Pest Spodoptera littoralis Caterpillars.

Starvation is frequently encountered by animals under fluctuating food conditions in nature, and response to it is vital for life span. Many studies have investigated the behavioral and physiological responses to starvation. In particular, starvation is known to induce changes in olfactory behaviors and olfactory sensitivity to food odorants, but the underlying mechanisms are not well understood. Here, we investigated the transcriptional changes induced by starvation in the chemosensory tissues of the caterpillar Spodoptera littoralis, using Illumina RNA sequencing. Gene expression profiling revealed 81 regulated transcripts associated with several biological processes, such as glucose metabolism, immune defense, response to stress, foraging activity, and olfaction. Focusing on the olfactory process, we observed changes in transcripts encoding proteins putatively involved in the peri-receptor events, namely, chemosensory proteins and odorant-degrading enzymes. Such modulation of their expression may drive fluctuations in the dynamics and the sensitivity of the olfactory receptor neuron response. In combination with the enhanced presynaptic activity mediated via the short neuropeptide F expressed during fasting periods, this could explain an enhanced olfactory detection process. Our observations suggest that a coordinated transcriptional response of peripheral chemosensory organs participates in the regulation of olfactory signal reception and olfactory-driven behaviors upon starvation.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization.

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

LRez: a C++ API and toolkit for analyzing and managing Linked-Reads data.

Motivation: Linked-Reads technologies combine both the high quality and low cost of short-reads sequencing and long-range information, through the use of barcodes tagging reads which originate from a common long DNA molecule. This technology has been employed in a broad range of applications including genome assembly, phasing and scaffolding, as well as structural variant calling. However, to date, no tool or API dedicated to the manipulation of Linked-Reads data exist. Results: We introduce LRez, a C++ API and toolkit that allows easy management of Linked-Reads data. LRez includes various functionalities, for computing numbers of common barcodes between genomic regions, extracting barcodes from BAM files, as well as indexing and querying BAM, FASTQ and gzipped FASTQ files to quickly fetch all reads or alignments containing a given barcode. LRez is compatible with a wide range of Linked-Reads sequencing technologies, and can thus be used in any tool or pipeline requiring barcode processing or indexing, in order to improve their performances. Availability and implementation: LRez is implemented in C++, supported on Unix-based platforms and available under AGPL-3.0 License at https://github.com/morispi/LRez, and as a bioconda module. Supplementary information: Supplementary data are available at Bioinformatics Advances online.